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From Feynman's path integral, we derive quasi-classical quantization rules in supersymmetric 
quantum mechanics (SUSY-QM). First, we derive a SUSY counterpart of Gutzwiller's formula, 
from which we obtain the quantization rule of Comtet, Bandrauk and Campbell when SUSY is 
good. When SUSY is broken, we arrive at a new quantization formula, which is found as good 
as and even sometime better than the WKB formula in evaluating energy spectra for certain one- 
dimensional bound state problems. The wave functions in the stationary phase approximation are 
also derived for SUSY and broken SUSY cases. Insofar as a broken SUSY case is concerned, there 
are strong indications that the new quasi-classical approximation formula always overestimates the 
energy eigenvalues while WKB always underestimates. 



I. INTRODUCTION 

The idea of supersymmetry (SUSY) is based on the 
expectation that there may be environments where the 
distinction between bosons and fermions is irrelevant 
[1]. Although there are some indications that the SUSY 
scheme works in understanding low energy phenomena 
[2], it is an observed fact that SUSY is generally broken. 
In 1981, Witten [3] utilized SUSY quantum mechanics 
to simulate spontaneous breaking of SUSY through non- 
perturbative quantum effects [4]. Despite its field the- 
oretic origin, SUSY quantum mechanics has many in- 
teresting properties that have been utilized for solving 
various problems in non-relativistic quantum mechanics 
[5]. In fact, as early as 1976, Nicolai [6] had already 
used the idea of SUSY quantum mechanics for the study 
of statistical mechanics. As a quantization method, it 
has been successful in reproducing exact solutions of the 
Schrodinger equation for a class of so-called shape invari- 
ant potentials [7]. This relatively new method has also 
proven useful in studying more complex problems such 
as tunneling problems [8], the Kaluza-Klein monopole 
[9], the three body partition function for an anyon gas 
[10], the Pauli and the Dirac Hamiltonian for electrons 
in magnetic fields [11,12] and the magnetic top [13]. 

In 1985, Comtet, Bandrauk and Campbell [14] pro- 
posed for SUSY quantum mechanics a semiclassical quan- 
tization formula, which henceforth we shall call the CBC 
formula, 
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where xl and xr are the left and right turning points 
given by (t>' 2 (x L ) = (t>' 2 (x R ) = E. This formul a is similar 



to but different from the well-known WKB quantization 
rule, 
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with V{x£) = V(xr) = E. In (1), E is the energy shifted 
so that E = for the ground state. It differs from E of (2) 
by the energy shift e, that is, E = E + e. The real- valued 
function <j){x) appearing in (1), which will be referred to 
as the superpotential by following the recent convention, 
is related to the potential V(x) of (2) by 

V(x) = <f> 2 (x)±-H=<f>'(x) + e, (3) 

where <j)'{x) = d(f)(x)/dx. Depending on the sign of the 
second term, E — V(x) equals E — V+(x) or E — V-{x) 
where V±(x) = (f) 2 (x) ± (h/\/2m) <j)'{x). The ground 
state is customarily assumed to belong to V-{x), so in 
(1) n = 0, 1,2, . . .for V-(x) and n = 1, 2, 3 . . . for V+(x). 

As is well-known, the WKB semiclassical quantization 
formula gives exact spectra for exactly soluble examples 
provided some Langer-like parameter modifications are 
made (see, e.g., [5]). Surprisingly, the CBC formula (1) 
has been found to be able to reproduce the same exact 
spectra without any ad hoc modification [14]. Even for 
those problems which are not exactly soluble, the CBC 
formula has been claimed to yield better approximation 
than the WKB formula does [15]. 

The success of the CBC formula is truly remarkable. 
Nevertheless, the origin of the formula is not very clear. 
There have been attempts to derive the CBC formula 
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from the WKB formula. Since the WKB formula itself 
has the Langer-like ambiguity, such attempts only verify 
that whenever the WKB formula is acceptable the CBC 
formula should work in the same semiclassical limit. It 
is certainly desirable to derive the CBC formula directly 
from what may be taken as the first principle. 

If SUSY is broken, however, the CBC formula has no 
ground for its validity. A natural question arises as to 
whether there will be a CBC counterpart for a broken 
SUSY case. In a previous work [16], we have proposed 
the following formula for a broken SUSY case, 



(4) 



with (/) 2 (xl) = <f) 2 {xR) = E. It has also been shown that 
this "broken SUSY formula" can reproduce exact spectra 
for the radial harmonic oscillator and the Poschl- Teller 
system [17,18]. 

In the present paper, we wish to show that the CBC 
formula (1) can be derived, when SUSY is good, from 
Feynman's path integral by a stationary phase approx- 
imation. Then we report that the CBC counterpart, 
applicable only to the case where SUSY is broken, can 
also be derived from the same path integral. Further- 
more, we propose that the new quantization formula for 
broken SUSY is of practical value as an approximation 
formula which is as good as or sometimes better than 
the WKB formula. However, we have no answer why 
the CBC formula can provide exact results for shape- 
invariant potentials and why the SUSY approximations 
whenever applicable are better than the WKB estima- 
tions. Section II briefly describes the backgrounds for 
our discussions. In Section III, we develop the quasi- 
classical approach to derive the quasi-classical counter- 
part of Gutzwiller's semiclassical formula for the energy- 
dependent Green function from which we obtain the CBC 
formula (1) and the broken SUSY formula (4). In Section 
IV, for a couple of examples with broken SUSY, including 
a class of power potentials, we analyze numerically the re- 
sults of the WKB formula and the broken SUSY formula 
in comparison with those from Schrodinger's equation. 
For the systems whose SUSY is broken, there are indi- 
cations that the standard WKB formula underestimates 
the energy values, while the broken SUSY formula over- 
estimates them. A way to improve approximation for the 
energy spectrum of a broken SUSY system is also sug- 
gested. Remarks are made in Section V concerning the 
limitation of the broken SUSY formula. 



II. SUSY QUANTUM MECHANICS AND PATH 
INTEGRALS 



Let us start with the SUSY invariant Lagrangian in 
one dimension [19]: 
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where t/> and t/>t are Grassmann variables. From this, 
Feynman's path integral can formally be constructed, 
which consists of integrations over the bosonic coordi- 
nate variables x as well as the fermionic Grassmann vari- 
ables t/> and t/>t . After carrying out path integration over 
the Grassmann variables [19], we can arrive at two path 
integrals in bosonic coordinate variables, 



x"=x(t) 



K± (x" , x'; t) 



J exp < — J L±(x,x)dt > T>x(t), 



x'=x(0) 



for the following effective Lagrangians, 
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L±(x, x) = -mi 2 — <f> 2 (x) =p 



-.<f>'(x). 



(6) 



(7) 



The last term in (7) is a non-perturbative quantum ef- 
fect arising from elimination of the fermionic degrees of 
freedom (the contribution of a//fermion loops). The La- 
grangians (7) correspond to the Hamiltonians, 



2m dx 2 



(8) 



which characterize the simplest version (N = 2) of SUSY 
quantum mechanics as defined by Witten [3]. Naturally, 
the path integrals (6) describe the propagators which 
are associated with the time-evolution generated by the 
Hamiltonians (8): 

K±(x", x'; t) = (x" | e - (i ' n > H ± | x'). (9) 
The Hamiltonians (8) may be factorized as 

H + = AA\ H- = A^A, (10) 

by means of the operators, 

k d , . . , + h 



A 



+ (j>{x), A^ 



+ <f>(x). (11) 



Therefore, it is evident that H± are positive semi- 
definite, i.e., spec(if±) > 0. For the eigenvalue prob- 



lems, H± ip^\x) = E^' ipl? 1 ' (x), we can show that 
E n > except for the ground state 



(±),/,(±V 



E 



(+) 



E 



(-) 



and that 



^+\x) = -JL^(-)(z), ^-\x) = -±= A^i+\x), 
V En V E n 



for E n ^ 0. Customarily, the ground state is as- 
sumed to be an eigenstate of H _ . Namely, by choice, 
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H- t/"o \ x ) = Eoipo \x). The lowest eigenvalue of H + 
is Ei. 

If SUSY is a good symmetry, that is, if the ground 
state is invariant under the SUSY transformation, then 



Aip[ \x) = 



(12) 



which means that there is no degeneracy in the ground 
state and that the ground state energy must vanish. 
Thus, when SUSY is a good symmetry, Eo = 0. This 
relation together with (11) also implies 



t/>o \x) = t/>o \x )exp 




(13) 



In fact, (14) is Feynman's path integral with a poten- 
tial of the form (16). However, the path integral (14), 
whether it is subjected to the condition (16) or not, is in 
general difficult to calculate. There are only a very lim- 
ited number of examples that can be exactly solved by 
Feynman's path integral. Historically, it is known that 
path integration for x 6 R can explicitly be carried out 
only for a quadratic system [20]. If we use polar coor- 
dinates, then we may include the harmonic oscillator in 
an inverse-square potential to the list of exactly soluble 
examples [21]. Furthermore, by mapping into a higher 
dimensional path integral, we can solve the Poschl- Teller 
oscillator and their variations [22,23]. It is interesting to 
note that many of these soluble potentials are reducible 
to the form (3). 



When SUSY is good, this solution must exist. In other 
words, the ground state function must be normalizable 
(or square-integrable). The normalizability of (13) poses 
a condition on the superpotential <f>(x). The ground state 
function (13) is normalizable only when the integral of 
the superpotential <j){x) tends to infinity as x goes to 
±oo. For instance, when we expect a purely discrete en- 
ergy spectrum, <j){x) must be such that |<^(a;)| —> oo as 
\x\ —> co. If <j){x) is a continuous function having an 
odd number of zeros and the ground state (13) has a 
vanishing eigenvalue, Eo = 0, then SUSY is a good sym- 
metry [3]. On the other hand, if <j){x) is continuous and 
has an even number of zeros, the corresponding ground 
state has a strictly positive eigenvalue and will be degen- 
erate. Hence SUSY is spontaneously broken. Therefore, 
the ground state of H- determines whether SUSY is bro- 
ken or not. For the present purposes, it is sufficient to 
consider _L_ of (7) which corresponds to H- of (8). 
Thus the path integral we wish to calculate is 



-x(t) 



K(x" , x' ; t) 



exp{(i/h)S[x(t)]}Vx(t) (14) 



=x(0) 



with the action for L. 



S[x(t)}= f (jx 2 -V.(x))dt, 



where 



V-(x) = (f> 2 (x) 



:<f>'(x). 



(15) 



(16) 



III. THE QUASI-CLASSICAL APPROACH 

For a more general potential system, we have to pur- 
sue an approximate solution. The most appropriate ap- 
proximation method in evaluating a path integral is the 
stationary phase calculation. In the usual semiclassi- 
cal stationary phase approximation the action functional 
^[^(t)] is expanded about the classical path x c \(t) as 

S[x(t)] ~ S d (x", x' , r) + 6S[x(t)] + 6 2 S[x(t)]. (17) 

The classical action S c \ = S c \(x" , x' , r) = S[x c \(t)] is an 
action evaluated along the classical path from x' to x" , 
determined by 8S[x(t)] = 0. Then the propagator is 
given by the formula of Van Vleck, Pauli and Morette 
(VMP) [24], 



K(x" , x' ; t) 



t fixed 

E 

x cl (t) 



2ttS dx' dx" 



exp 



-r^ci 
n 



(18) 



If there are classical paths more than one, we have to 
add up the contributions from all these paths. Thus, 
the summation in (18) is to cover all the classical paths 
x c \(t) between x' and x" with a fixed time interval r. 
This approximation formula is known to give rise to the 
exact propagator for the free particle, the harmonic oscil- 
lator and more general quadratic systems, but does not 
directly provide the WKB quantization rule (2). The 
energy-dependent Green function evaluated by summing 
over the classical paths for a fixed energy E is the source 
of the WKB formula. The Fourier transformation of the 
semiclassical propagator (18) leads to Gutzwiller's for- 
mula [25], 



E fixed 



G(x",x';E) ~ — 



„.. J2 exp{(,yft)W cl -ifr/2M ; 



(19) 
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where W c \ = W c \ (x" , x' , E) = S c i + Et is Hamilton's 
characteristic function evaluated along the classical paths 
x c \(t) with fixed energy E, and v is an integral num- 
ber known as the Morse index. The WKB quantization 
rule (2) results from the poles of this Green function. 
Again, the WKB formula derived by this semiclassical 
approximation gives the exact energy spectrum for the 
harmonic oscillator in one dimension. However, it does 
not yield the exact result for the three-dimensional ra- 
dial harmonic oscillator unless the Langer replacement 
1(1 +1) -> (I + 1/2) 2 is made. For other exactly 
soluble examples, exact spectra may be obtained from 
the WKB formula by applying appropriate Langer-like 
modifications. In other words, most WKB results for ex- 
actly soluble systems are not exact without Langer-like 
replacements. The Langer replacement and other simi- 
lar modifications are ad hoc procedures which reflect the 
ambiguity involved in the choice of terms of 0(h 2 ) for 
a stationary action. Applying a Langer-like replacement 
amounts to introducing a so-called "quantum correction" 
term V c (x) in the Lagrangian. This means that an ap- 
propriate effective action, 



S efl [x(t)] 



V.(x)-V c (x)) 



dt, 



may be chosen so as to generate a result corrected by 
a Langer-like modification. So far, no general rules 
have been established for the choice of such correction 
terms. Therefore, whenever a Langer-like modification is 
needed, there is no inherent reason why we should make 
the action of the form (15) stationary. In fact, we are 
free in principle to make any action stationary. Since 
our desire is to improve the semiclassical approximation 
within the framework of the stationary phase approach, 
we must look for a more suitable effective action which 
may be made stationary. 



A. Sum over Quasi-Classical Paths 

In order to seek an alternative and equally qualified 
approach, let us recall that the derivative term in (16) 
is the contribution from all fermion loops, while the re- 
maining part of the action (15) contains only the contri- 
butions of the tree diagrams in perturbation calculations. 
It is uncertain how much such background information 
will remain relevant when we accept the SUSY quantum 
mechanical method as a tool in non-relativistic quantum 
mechanics. Nonetheless, we consider that the key to the 



alternative approach lies in separating the action (15) 
into two parts, 



S[x(t)] = S tTee [x(t)] + h<p[x(t)], 
with the "tree" action, 



(20) 



Street)] = / (j ~ ^ ( X )) dt > ( 2 1) 



and the "fermionic" action representing the fermion loop 
corrections, 



<p[x(t)] 




<f>'(x(t))dt. 



(22) 



Then, we pursue an approximation for which the tree 
action (21) is made stationary. 
Namely, we demand that 

aStree [?(*)] = ( 2 3) 

which determines specific paths x = q(t) with momentum 

p qc (q) = mq = ±^<2m(E - ^(q)), (24) 

and energy E = — dSt lTee [q(i)\l dr . Here, the + and — 
signs are to be taken for the motion to the right and left, 
respectively. Since the "classical" paths x c \(i) are speci- 
fied by 8S[x(t)] = 0, the paths q(t) along which the tree 
action (21) is stationary are not quite classical. We shall 
call them "quasi-classical" paths, and we shall carry out 
the quasi-classical stationary phase calculation 1 for the 
path integral (14). 

Now we expand the tree action (21) about a quasi- 
classical path q(t) to second order by letting £(t) = 
x(i) — q(i) with x(0) = q(0) = a? c i(0) = x' and x(t) = 
q(r) = x c1 (t) = x": 



S tlee [x(t)} ~ S qc (x", x', t) + 6 2 S[£(t)] 



(25) 



qc = S qc (x" , x' , r) = Stree[q(t)] is the tree action 



Here, S ( 

(21) evaluated along the quasi-classical path q(t), and the 
second variation 8 2 S[£,(t)], representing the contribution 
from the fluctuation about the quasi-classical paths, is 
given by 



— rat 

2 ^ 



2 T 2 mn 2 (t)e 



dt 



(26) 



where mQ?{t) = d 2 <f> 2 (q(t)) / dq 



1 In the literature the terms "semiclassical" and "quasi-classical" are often synonymous. Here, in contrast, we are differentiating 
the quasi-classical calculation from the usual semiclassical approximation. 
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Then we calculate the path integral (14) by the station- 
ary phase approximation. Because of the quasi-classical 
path dominance, the fermionic action (22) takes the form, 



f = <pW)} 



4>( X ") 



±\/E- 



(27) 



</>(V) 



which immediately integrates, if no turning points are 
involved, to 



<p = \ (a(x") - a(x')) 



where 



a(q) 



(28) 



(29) 



In general, the quasi-classical motion takes place peri- 
odically between the left turning point xl and the right 
turning point xr of the quasi-classical paths. At these 
turning points, p qc {xL) = Pqc{%R) = 0, and the momen- 
tum changes its sign. The fermionic action (22) thus 
turns to an overall phase to the partial propagator cor- 
responding to each quasi-classical path: 



<Pk = <Po + k[a(x R ) - a(x L )] 



(30) 



where k signifies the number of full cycles of the path 
q(t), and ipo is the remaining phase due to incomplete 
cycles. If, in particular, the path has no turning point, 
the phase takes the special form (28). 

The quasi-classical approximation of the action (20) 

S[x(t)} ~ 5 qc + hip + 6 2 S[((t)} 

reduces the path integral (14) to 

t fixed 



K(x" , x' ; t) 



J2 A'(0,0;r)expj-5 < 
«(*) 



qc 



lip 



(31) 



where 



e(r) = 



#(0,0; r) 



exp { -8 2 S[C(t)] } Wit) ■ 



(32) 



e(o)=o 



The remaining path integral (32) for the second variation 
(26) over the ^-variable can easily be calculated by the 
standard technique [26], which results in the prefactor of 
the form, 



#(0,0; r) 



2™/i/("r) 



(33) 



The function f(t) appearing in this prefactor is the solu- 
tion of 



dt 2 



f(t) = 0, 



(34) 



which satisfies the initial conditions, /(0) = and /(0) = 
1. The solution can be put, as is easily checked, in the 
form (see, e.g., [27]), 



fir) = 3(r),j(0) 



dt 



WW 



(35) 



Since the velocity q(t) changes sign at every turning 
point, we have the relation, 



m 2 q(T)q(0) = (-iy\p qc (x")p qc (x')\ 



(36) 



Here v is the number of turning points along a particular 
quasi-classical path q(t) between x' and x" , that is, the 
number of the zeros of f(t) for t £ [0, r], which is in this 
case the same as the Morse index [28]. It is zero when r 
is small, and hence the solution /(t) of (34) satisfies the 
desired initial conditions. 

The right hand side of (35) may also be expressed in 
different fashions. Integrating the inverse of (24) gives 
the time interval r as a function of d E, 



(37) 



±^2m(E-e(q)) 



In the last equality, the negative sign is to be chosen when 
q < 0. Differentiating both sides of this yields 



dE 



1 f dt 

m 



o 



WW 



(38) 



Combining (36) and (38), we may rewrite (35) as 



f(T)=--(-iy \p q c(x")p qc (x 

m 



dE 



(39) 



If we further define the quasi-classical counterpart of 
Hamilton's characteristic function, 



W qc (x",x',E) = S qc (x",x',T) + ET= / Pqc (q) 



we get 



dW, 



dE 



qc d 2 W qc 3t 

T, - 



dE 2 dE 



(40) 



(41) 



as well as 



Pqc(x") 

Pqc(x') ~- 



dW qc _ 


dS qc 


dx" 


dx" ' 


dW qc 


_ d,S qc 


dx' 


dx' 



(42) 
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At this point, if we wish, we can derive the quasi- 
classical counterpart of the VMP formula (18) in a way 
parallel to the standard derivation of (18). From (42) it 
is obvious that 



fixed r 



«(*) V 



d 2 S< 



qc 



dx' dx" 



exp 



-5 qc + if - i—v 



(43) 



d 2 S qc _ d Pqc (x") 



dE 



dx' dx" 



dx' 



Pqc(«") 9x' 



where we have used the relation dp qc (q)/dE = m/p qc (q) 
resulting from (24). From the Hamilton-Jacobi equation 
for the quasi-classical paths, E = —d,S qc /dT, there fol- 
lows 



dE 
dx' 



d 2 S qc dp qc (x') m 



dr 



B. Gutzwiller's Formula for SUSY-QM 

To find the quasi-classical counterpart of Gutzwiller's 
formula (19), we have to go over from the propagator 
to the energy-dependent Green function by the Fourier 
transformation, 



drdx' dr Pqc(x') \dE t 

Using these results, we replace dr/dE in (38) by 



G(x", x'; E) = — / K(x", x'; r) e ( ^ n ^ T dr. (44) 

In going over to Gutzwiller's semiclassical formula (19), 
the VPM formula (18) is usually used [26,29]. In much 
the same manner, we may utilize the quasi-classical coun- 
terpart (43) of the VPM formula to compute the time- 
and arrive at a quasi-classical counterpart of the VPM integral of (44). However, for clarity, we choose a slightly 



dr\ 1 1 , „ 5 2 5 qc 



dx' dx" 



formula, 



different approach. 



Substituting (39) and (41) into (33), we express the propagator (31) in the form, 

t fixed I 



K(x",x';t)~ J2 



1 



«(*) 



27riR(-l)"|p qc (a:')Pqc(a;")l 



d 2 W a , 



dE 2 



-1/2 



exp{-[W qc (E)- Er] + if 



(45) 



Then we put this into (44) and convert the r-integration into the £"-integration by letting dr = \dr/dE \ dE 
\d 2 W qc /dE 2 \dE to get 



G(x",x';E) 



E 



hj ^ V 2irih(-iy\p qc (x')p qc (x")\ 

1{ t ) 



d 2 W a , 



dE 2 



1/2 



exp { -W qc (E) + -(E-E)t + if } dE. (46) 



Next, we expand the function W qc (E) about E as 



Wqc (E) = W qc (E)+(^ 



( F F )+ l f d2W ^ 



(E-E) 2 



(47) 



and perform the stationary phase integration over E about E. In this manner, we are able to arrive at the following 
quasi-classical approximation formula for the energy-dependent Green function, 



G(x" , x'; E) 



1 



E fixed 



ih l \ ( '\ i //M ^ exp \ii/Ti)W qc {E) + if - i(7r/2)^} 



(48) 



where the sum is over all quasi-classical paths q(t) from x' to x" with a fixed energy E. This is the quasi-classical 
counterpart of Gutzwiller's formula (19). 
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C. Quasi-Classical Quantization 



w(x) = I ^2m(E - <f> 2 (q))dq, (50) 



XL 



For each quasi-classical trajectory with the same fixed 
energy value E, we have to find W qc , <p and v in order 

to perform the summation in (48). Following Schulman's and a ^ =a , IC smU(x)/\^E) as given in (29). The quan- 



prescription [261, we group the set of all paths into four , „ , P > 

, L J titles defined for k = are 

classes: 

(1) paths which leave x' to the right and arrive at x" Wq = w(x") — w(x'), 

fr0mtheMt; <pW = l[ a (x")-a(x')], 

(2) those which leave x' to the right and arrive at x" ^(l) _ q 
from the right; ' 

(3) those which leave x' to the left and arrive at x" 

from the left; and i/^ 2 ) = w ( x ") + w(x'), 

(4) those which leave x' to the left and arrive at x" = i [a(x") + a(x')] — a(«i) 
from the right. , . 

v a = l > 

They may be reflected at the left turning point xl, or 
the right turning point x R , or at the left and right turn- 
ing points. Within each class, denoted by the superscript w( 3 ) 
i (i = 1, 2, 3, 4) in parentheses, a path is uniquely charac- 
terized by a non-negative integer k indicating the number ^(3) _ a(x R ) — - [a(x") — a(x')] 
of full cycles the path contains. For a path with k cycles 2 
in class (i) (i = 1,2,3,4), we write (W^ ) = W ( k i) ) ^ = 1 , 



W [ o> = 2w(x R )-w(x")-w(x'), 



W ( k i] = + 2kw(x R ), and 

= <Po } + k H x n) ~ «<>£)] > ( 49 ) = 2w(x R )-w(x") + w(x'), 



} = v o } + 2*, ^ = a (x R ) - a(x L ) - \ [a(x") - a(x')] 

where j^ 4 - 1 = 2. 



Now the sum for a fixed E in (48) can be rewritten as 



E fixed 4 oo 

E (•) = £»•)• 

q(t) i = l k = 



The second summation leads to a geometric series which can be easily evaluated. As a result, the Green function in 
the quasi-classical approximation becomes 



E exp {(i/n)W^ + - i(T/2)^°} 



G(x" x'- E) ~ . i=1 (51) 

ihy/\ Pqc (x') Pqc (x")\ 1 - exp {i [2w(x R )/h + a(x R ) - a(x L ) - tt]} ' 

Apparently, the poles of the Green function (51) occur when 

w(x R ) = in + — 1 jrh (52) 

where n = 0,1,2,.... This is indeed a quasi-classical quantization condition, which implies two different formulas 
depending on whether SUSY is good or broken. To determine the values of a(xj^) and a(x R ) explicitly, we have to 

recall the conditions at the two turning points, <f> 2 (xL) = <f> 2 (x R ) = E, or |«^(a? _t) | = |<^(#ii)| = which has two 

solutions [16]. Accordingly, we have the following two distinct cases: 
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Case I: 4>(xr) = — 4>(xl) = ivE. In this case, it is apparent from (24) that ci(xr) = —a(xL) = ± arcsin 1 = ±7r/2. 
Since the ground state is assumed to be an eigenstate of H-, the upper sign should be selected. 2 The fermionic phase 
ip has a non-zero contribution which cancels the fractional part on the right hand side of (52). Consequently, the 
quantization condition (52) results in 

w(xr) = nirh, n = 0,l,2,..., (53) 

which coincides with the celebrated formula (1) of Comtet, Bandrauk and Campbell [14]. Since <j)[x) must have an 
odd number of zeros in this case, SUSY is a good symmetry. For the residues of the Green function G(x" , x' ; E) at 
the poles due to (53), we have 



ResG(x",x';E) 

where 



E=E n 



4m cos (j- w(x') + \ a(x')) cos (j- w(x") + \ a(x")) 



2m/(E n -^(q))dq (55) 



is the period (the time interval for a single cycle) of the bounded motion with energy E n . From this follow the excited 
state wave functions for Ql < % < 1r, 

, /4m , , , , -i /, / w(x) 1 <t>(x) \ 

^\x)~ J— p c i a: " 7 cos ^Ttt arcsin ) , (56) 

which is normalized within the quasi-classical scheme. In the path integral approach, the results are usually normalized 
by the basic requirement, 

lim K(x" , x'; r) = 6(x" — x'). 

Note that these wave functions are only for the excited states with E n > (n = 1, 2, 3, . . .). The ground state wave 
function has already been given by (13). 



Case II: 4>{xr) = <f>{%L) = ±V E. In this case, cl{xr) = cl{xl) = ± arcsin 1 = ±7r/2 and the contribution of the 
fermionic phase <p vanishes in (52). Therefore, the quantization condition (52) implies another formula, 

w(x R ) = (^n + irh , n = 0,l,2,..., (57) 

which is identical with the formula (4) we have proposed for the broken SUSY case. The left-hand side of this formula 
is the same as that of the CBC formula (1), whereas the right-hand side equals that of the WKB formula (2). This 
new formula, which is a hybrid of the WKB and the CBC formula, is nothing but the usual WKB result for the 
"tree" Hamiltonian without the fermion loop contributions (ip = 0). Since <j)[x) has an even number of zeros, SUSY is 
broken. The hybrid formula (57) is indeed the CBC counterpart for broken SUSY cases. Although the fermion loops 
do not affect the spectrum in the present approximation, there is an effect on the quasi-classical eigenfunctions. The 
residues of the Green function at the poles due to (57) are 



Res G(x" , x' ; E) 



4m sin (i w(x') + \ a(x')) sin (± w(x") + \ a(x")}_ 



E=E n T n \J\Vc\{x')\ VM*")I 

We immediately find the quasi-classical wave functions, 



2 To obtain the quasi-classical spectrum for H+ the lower sign has to be chosen. This leads to the same quantization condition 
(53). However, n = 1, 2, 3, . . . for this case. 
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T I arcsm iM] ; (59) 
2 \/E n 



where n = 0,1,2,..., and < x < qr. 



The quasi-classical quantization condition (52) derived in the above from the path integral is identical with the 
rule obtained earlier by Eckhardt [31] from Maslov's asymptotic analysis. In fact, he also identified the first case (53) 
with the CBC formula. However, he did not recognized that the second case (57) corresponds to broken SUSY. The 
hybrid formula (57) was for the first time related to broken SUSY in ref. [16] and applied to solve quantum mechanical 
problems in ref. [17,18]. 



IV. APPLICATIONS OF THE BROKEN SUSY 
FORMULA 

A. Exactly Soluble Examples 

It is a remarkable fact that the CBC quantization rule 
(1) supplies exact energy spectra for shape invariant po- 
tentials while the WKB rule (2) needs the Langer-like 
modification. What may as well be remarkable is that 
the broken SUSY formula (57) can also provide without 
any ad hoc modification the exact spectra of systems for 
which the superpotential has a general form [17,18], 

<f > (x) = Af(x) + B/f(x), (60) 

where A and B are constants. The radial harmonic os- 
cillator, the Poschl-Teller oscillator, and the modified 
Poschl-Teller oscillator are such examples [17,18]. Since 
these examples are basic, it would be instructive to com- 
pare SUSY compatible forms and broken SUSY compat- 
ible forms of the superpotential explicitly. 

(a) Tie radial harmonic oscillator: 

V.(r) = \muj\ 2 + /(/ + y + (r, - l/2)hu , 
2 2mr l 

(1 = 0,1,2,..., r> 0) for which 

[m Ti 

n x ) = \hr ur + V — , 
V ^ v 2m r 

where rj = — / — 1 (good SUSY) or rj = / (broken 
SUSY). 

(b) The Poschl-Teller oscillator: 

V- (x) = Vo [k(k — 1) esc 2 x + A(A — 1) sec 2 a;] 
-V ( V -Xf, 

(k > 1, A > 1, < x < 7r/2) for which 

4>(x) = \/Vo [A tan x + rj cot x] , 

where rj = — k (good SUSY) or rj = k — 1 (broken 
SUSY). 



(c) The modified Poschl-Teller oscillator: 

V-(x) = Vo[k(k — 1) csch 2 a; — A(A — 1) sech 2 a;] 
-V ( V -Xf, 

(k > I, X > I, x > 0) for which 

4>{x) = [>tanh(^) + r]coth(x)], 

\/2m 

where rj = — k (good SUSY) or rj = k — 1 (broken 
SUSY). 

There are other systems soluble exactly by the broken 
SUSY formula [32], but those soluble are reducible to 
the form given above. 



B. The Power Potential Problems 

The new quasi-classical formula (4), like the WKB for- 
mula (2) and the CBC formula (1), can serve significantly 
as an approximation formula even for non-shape invari- 
ant systems. In order to present a conspicuous difference 
between the approximate results of a SUSY system and 
a broken-SUSY system, let us studies a few examples. 

The first example we consider is a class of superpoten- 
tials of the form, 

<j,(x) = Ax k , A>0, k = 1,2,3,... , (61) 

which is not shape-invariant and not exactly solvable. In 
a recent work [15], Khare has studied the case of the odd 
power, that is, the case of A; = 2d — I where d is a positive 
integer. Since the superpotential is of an odd power, it 
has an a single zero but a zero of order 2c? — 1. There- 
fore, the CBC formula is applicable. In the case of an 
even power, say, k = 2d, the zero of the superpotential 
is of order 2d, so that SUSY is broken. The CBC for- 
mula is no longer applicable. In the earlier work [15], 
the even power case has been left out from being treated 
within the framework of SUSY quantum mechanics. The 
new broken SUSY formula (4) enables us to handle even 
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power cases and to supplement Khare's analysis of the 
odd power problem. 3 

The partner potentials derived from (61) for k = 2d 
are 



V±(x) = A 2 x 4d ± 



2hAi 



'2d-l 



(62) 



Obviously, V+{x) = V-(—x). Hence it is sufficient to 
consider V(x) = V-{x) + e. As is the case of the CBC 
integral for an odd power [15], the integral in the new 
broken SUSY approximation formula (4) can be explic- 
itly calculated, yielding an approximate energy spectrum 
expressed in a closed form, 



E n 



A 2/(2d+l) ( 

2m 



2 x 2d/(2d+l) 



4d/(2d+l) 



(63) 



We compare numerically these energy values with 
those obtained from the usual WKB approximation (2). 
In addition, we make exact numerical integrations of 
the Schrodinger equation. The relative deviations A = 
(-^exact - -E'approxV-E'exact (in %) from the exact eigenval- 
ues are shown graphically in Figure 1 for d = 1, 2 and 3. 
We have chosen units such that A = m = Ti = 1. For the 
ground state energy shift we have taken e = 2. Graphs 
of the corresponding potentials V_ (x) are given in Figure 
2. 

The results for d = 1 show that the WKB approxima- 
tion of the ground-state energy (n = 0) is better than 
the estimate by formula (4). Even for excited states the 
results of WKB are generally better than those from (4). 
However, for the cases where d = 2 and 3, the broken 
SUSY formula (4) provides better approximations for the 
ground-state energies. Except for the first few exited 
states, formula (4) is always better than WKB. In the 
limit d — ► oo, we obtain from (63) 



lim E n = (h 2 Tr 2 /32m) (2n 

d— ^oo 



i) 2 



(64) 



which coincides with the energy spectrum for the even 
states of a particle in an infinite square well of width 4 
(— 1 < x < 3). In fact, for large d, the potential V-{x) 
has the form of an infinite square well of width 2 with 
an additional deep negative dip near x = 1. The area 
enclosed by the negative part of V_ (x) and the positive 
a;-axis diverges as d — ► oo. This is in contrast to the odd 
power potential considered by Khare, which reduces in 
the case of V+(x) to a standard one-dimensional box in 
the limit d — ► oo. 



The energy eigenfunctions for large d, as observed by 
numerical integrations, are of the form 



ip n (x) = N F(x) sin [(2n + l)(x + 1)tt/4] 



(65) 



where N is a normalization constant and the envelope 
function F(x) has for large d a shape as shown in Fig- 
ure 3. In the limit d — ► oo it consists of two unit-step 
functions denoted by 6: 



lim F(x) = \ (e(x + 1) - 9{x - 1)) . 



(66) 



This means that the potential V_ (x) becomes that of an 
infinite square well with Dirichlet boundary conditions 
at x = — 1 and Neumann boundary conditions at x = 1. 
It should be noted that the behavior of V_ (x) near x = 1 
is very similar to that of the derivative-delta potential 
U (x) = —f3b'(x — 1) where [3 — ► oo as d — ► oo. It is known 
that U (x) for [3 — ► oo produces the Neumann condition 
at x = 1 and decouples the two half-lines x < 1 and 
x > 1 [33]. These asymmetric boundary conditions are 
in a sense remnants of the broken SUSY. Therefore, in 
the limit d — ► oo, the potential V(x) becomes equivalent 
to the left half (—1 < x < 1) of an infinite square well for 
— 1 < x < 3 having the (^'-potential at x = 1, in which 
only the even states of the standard infinite square well 
of width 4 survive. The spectrum (62) is indeed exact 
in the limit d — ► oo. This would suggest that for large 
d the new formula (4) will give better energy estimates 
than the WKB formula (2) does. 

Although the broken SUSY formula (4) does not nec- 
essarily improve approximations as is seen in the case of 
d = 1, it adds an important advantage if applied in com- 
bination with the WKB formula (2). A closer look at 
the relative errors given in Figure 1 makes us realize that 
the WKB estimates are always below the exact values 
and that the results of the broken SUSY formula always 
above the exact values. The mean values of the results 
of (2) and (4) will give improved energy estimation. This 
interesting feature is not unique to the power potentials. 
There are strong indications that for any broken SUSY 
system the WKB result is an underestimation whereas 
the broken SUSY calculation results in an overestima- 
tion. 



C. More Examples and Improved Approximations 

Let us examine additional examples for which the su- 
perpotentials are not power functions: 



(a) <j){x) = cosh x and 



3 After the first version of this analysis was submitted for publication, a similar work on the even power potential appeared in 
ref. [30]. The emphasis there differs from ours. 
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the examples (a) and (b). In both cases the broken SUSY 
formula is better than the WKB formula. This is clearly 
demonstrated by the negative sign before the deviations 
of the improved values from the exact results. Hence 
averaging of the WKB and the CBC results has no par- 
ticular significance. In Figure 5, we show a numerical 
comparison of the WKB estimation and the CBC approx- 
imationfor the odd potentials <f>(x) = x\x\, <f>(x) = x 3 and 
<f>(x) = x 5 . This time, we employ the shift e = 2. Natu- 
rally, the WKB formula (2) yields different results for the 
two potentials V-{x) and V+{x). In Figure 5, although 
we use the same symbols for both WKB results, the data 
for V+{x) are connected by thick lines while those for 
V-(x) are connected by thin lines. For the CBC formula 
(1) or (52), we have to use by choice only V_ (x) to which 
the ground state belongs. Note that the lowest energy 
for V+(x) is Ei. Since the ground state energy of V-{x) 
vanishes in this case, the relative deviations A become 
infinite. Therefore, the estimates for the ground state 
energy by the WKB formula are omitted in Figure 5. 
As is apparent in Figure 5, the WKB calculation for V_ 
sometimes overestimates and the CBC value may become 
lower than the Schrodinger evaluation. The CBC approx- 
imation for <j){x) = x\x\ with n = 2 gives A = —0.009% 
and thus is slightly below the exact value. 



TABLE I. Comparison for (f>(x) = cosh x. The relative errors parenthesized are given in % . 



n 


exact 


WKB 


broken SUSY 


improved values a 





1.7190 


1.6688 


1.7650 


1.7169 






(-2.92) 


( + 2.68) 


(-0.12) 


1 


3.5578 


3.5165 


3.5890 


3.5528 






(-1.16) 


( + 0.88) 


(-0.14) 


2 


5.7122 


5.6754 


5.7380 


5.7067 






(-0.64) 


( + 0.45) 


(-0.10) 


3 


8.1421 


8.1088 


8.1652 


8.1373 






(-0.41) 


( + 0.28) 


(-0.06) 


4 


10.8210 


10.7906 


10.8430 


10.8168 






(-0.28) 


( + 0.20) 


(-0.04) 


5 


13.7291 


13.7021 


13.7516 


13.7269 






(-0.20) 


( + 0.16) 


(-0.02) 


a Obtained by 


averaging the WKB and broken SUSY result 


according to (67). 






TABLE II. Comparison for (f>(x) 


= exp(x 2 /2 


). The relative errors parenthesized are 


given in % . 


n 


exact 


WKB 


broken SUSY 


improved values 





1.7663 


1.6684 


1.7965 


1.7325 






(-5.54) 


( + 1.71) 


(-1.92) 


1 


3.8391 


3.7455 


3.8712 


3.8084 






(-2.44) 


( + 0.83) 


(-0.80) 


2 


6.4793 


6.3874 


6.5126 


6.4500 






(-1.42) 


( + 0.51) 


(-0.45) 


3 


9.6317 


9.5414 


9.6664 


9.6039 






(-0.94) 


( + 0.36) 


(-0.29) 


4 


13.2606 


13.1727 


13.2976 


13.2352 






(-0.66) 


( + 0.28) 


(-0.19) 


5 


17.3411 


17.2562 


17.3810 


17.3186 






(-0.49) 


( + 0.23) 


(-0.13) 



(b) <f>(x) = exp(a; 2 /2). 

Both of these superpotentials satisfy the symmetric con- 
ditions <f>(xL) = 4>(xr) = \I~E~ at the turning points. 
Naturally, the CBC formula is inapplicable. Thus, we 
calculate the energy values by using the WKB formula 
and the broken SUSY formula, and compare their results. 
Here, for convenience, the ground state energy shift has 
been set to e = 0, which implies E = E. In Figure 
4, we observe again that the WKB formula (2) under- 
estimates the energy values whereas the broken SUSY 
formula (4) overestimates them. We have also examined 
other examples [18,34]. Remarkably, all the examples we 
have examined so far have commonly exhibited the fea- 
tures that the results of the broken SUSY formula are 
as good as the WKB approximations and that the two 
estimations are separately distributed above and below 
the exact values. Therefore, averaging the results from 
the two formulas will lead to an improved approximation 
for a system whose superpotential satisfies the condition 
(j>(x L ) = (f)(x R ) = ±\/~E: 

Eimp = \ (-EwKB + -^broken) • (67) 

In Table I and Table II, numerical comparisons between 
the improved values and the exact values are shown for 
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D. A Nonlinear Field Model 

Let us also suggest a possible application of the bro- 
ken SUSY in field theories. In analogy to the class of 
potentials characterized by the superpotential (62) we 
may consider a one-dimensional field theory defined by a 
Lagrangian density of the form 



dx 



/?$ 4d + 7tft— + n$ 2d - 1 [tf t ,tf] (68) 
dx 



which characterizes a one-dimensional Ginzburg-Landau 
field <&(x) coupled to a fermionic field ^(x). It is known 
[36] that the inverse correlation lengths for such models 
are given by the energy eigenvalues of the corresponding 
Hamiltonian (8). A quasi-classical estimate of the eigen- 
values by (2) and (4) will thus give simple lower and 
upper bounds of correlation lengths, respectively. This 
might lead to some insights, in particular, to the nature 
of phase transitions of such field theories even in more 
than one dimension. A recent study of Kapuscik, Uzes 
and Barut [37] shows that an exact result of the nonlin- 
ear complex oscillator yields a spectrum very similar to 
the present quasi-classical result. 



V. CONCLUDING REMARKS 

In this paper, we have derived from Feynman's path 
integral the SUSY counterparts of the Van Vleck-Pauli- 
Morette formula for the propagator and the Gutzwiller 
formula for the energy dependent Green function by the 
quasi-classical approximation. From the poles of the 
Green function in the quasi-classical limit, we have ob- 
tained the CBC formula (1) for exact SUSY and the new 
quantization condition (4) for broken SUSY. We have 
then pointed out a limited role of the broken SUSY for- 
mula in reproducing exact energy spectra and demon- 
strated by few examples that the broken SUSY formula, 
whenever applicable, is generally as good as and some- 
time better than the WKB formula. We have also ob- 
served that the interesting feature that in approximation 
the broken SUSY formula overestimates while the WKB 
formula underestimates. Because of this special feature, 
the approximation can be improved by taking an average 
of the results from the two formulas. 

Unfortunately, the present approach stops short in an- 
swering the question as to why the CBC formula, and 
the broken SUSY formula in a limited extent, can yield 
exact spectra for shape-invariant potentials without the 
Langer-like modification. It has been argued that the 
exactness of the lowest stationary phase approximation 
may be assured if all the higher order corrections are 
zero. However, it is difficult to show in a general way 
that all the higher order contributions vanish. Further- 
more, there are instances that the first few order cor- 
rections even to the WKB approximation vanish while 



the Langer-like replacement is still required. It has never 
been proven that the Langer-like correction to the WKB 
result is indeed stemming from the higher order terms in 
k. The question remains to be answered. 

As has been discussed in the end of Section III, the con- 
ditions at the turning points helps to discriminate good 
SUSY and broken SUSY. The antisymmetric condition 
4>(xl) = —4>{xr) corresponds to a case where SUSY is 
good, while the symmetric condition 4>{xl) = 4>{ x r) lm ~ 
plies broken SUSY. If <f) 2 {x) is symmetric about the axis 
of a; = (xl + xr)/2, then the parity of the superpoten- 
tial immediately dictates whether SUSY is broken or not 
[16]. It may be worth noticing that the fermionic phase 
ip for full cycles precisely cancels the Morse index for 



each path, that is, ip^ — ip^ = (tt/2) (v^ 

SUSY is good. This is not true when SUSY is broken. 
Recently, it has also been pointed out that there is a 
close connection between the Maslov corrections and the 
parity of the underlying problem [35]. 

Finally, we have to touch upon shortcomings of the 
quasi-classical approach. For a shape-invariant poten- 
tial, the scalar function <j){x) can be found in a rather 
simple form. However, it is not easy to find the function 
4>{x) for an arbitrary potential V(x). Even if we can get 
the function in an analytical form, the integration of the 
CBC formula (1) or the new formula (4) will be diffi- 
cult. For instance, if V(x) = ax", we have the Riccati 
equation [38] to solve 



,(0 _ „(0 



when 



(f>'(x) + a<f> 2 (x) = bx", 



(69) 



with a = —\/ c 2m/Ti and b = —ay/2m/h, whose solution 
is given by <f>(x) = a~ x u' \x) / 'u(x) with 



u(x) = y/x Z l 



2iv ab 



,(n + 2)/2 



(70) 



where Z v (z) is a solution of Bessel's differential equation; 
for example, the Bessel function J v (z). Therefore, quan- 
tization cannot be done analytically by either the CBC 
formula or the broken SUSY formula (4). The WKB for- 
mula (2), however, leads to analytical results in this par- 
ticular example. For a more general potential, we must 
solve the generalized Riccati equation, 



(f>'(x) + a<f> 2 (x) = V(x), 



(71) 



which is not an easy task. In this regard, neither the 
CBC nor broken SUSY formula surpasses the breadth 
of the WKB formula as an approximation formula even 
though there are occasions where the SUSY formulas give 
better estimations. Yet, we have to recognize that for 
many shape-invariant potentials the SUSY formulas have 
a mysterious power to yield exact results which the WKB 
rule fails in providing. 
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FIGURES AND FIGURE CAPTIONS 
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FIG. 2. The potential V-(x) in (62) for parameters d = 1, 2 and 3. 
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